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Abstract 

In  this  paper,  yield  functions  describing  the  anisotropic  behavior  of  textured  metals  are  proposed. 
These  yield  functions  are  extensions  to  orthotropy  of  the  isotropic  yield  function  proposed  by 
Cazacu  et  al.  (Cazacu,  O.,  Plunkett,  B.,  Barlat,  F.,  2006.  Orthotropic  yield  criterion  for  hexagonal 
close  packed  metals.  Int.  J.  Plasticity  22,  1171-1194).  Anisotropy  is  introduced  using  linear  transfor¬ 
mations  of  the  stress  deviator.  It  is  shown  that  the  proposed  anisotropic  yield  functions  represent 
with  great  accuracy  both  the  tensile  and  compressive  anisotropy  in  yield  stresses  and  /--values  of 
materials  with  hep  crystal  structure  and  of  metal  sheets  with  cubic  crystal  structure.  Furthermore, 
it  is  demonstrated  that  the  proposed  formulations  can  describe  very  accurately  the  anisotropic 
behavior  of  metal  sheets  whose  tensile  and  compressive  stresses  are  equal. 

It  was  shown  that  the  accuracy  in  the  description  of  the  details  of  the  flow  and  /--values  anisotropy 
in  both  tension  and  compression  can  be  further  increased  if  more  than  two  linear  transformations 
are  included  in  the  formulation.  If  the  in-plane  anisotropy  of  the  sheet  in  tension  and  compression 
is  not  very  strong,  the  yield  criterion  CPB06ex2  provides  a  very  good  description  of  the  main  trends. 
©  2007  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Characterization  and  modeling  of  the  anisotropy  in  the  plastic  response  of  metals  with 
cubic  crystal  structure  is  well  advanced.  Anisotropic  yield  functions  that  capture  with 
increased  accuracy  both  the  anisotropy  of  the  yield  strengths  and  the  anisotropic  distribu¬ 
tions  of  the  Lankford  coefficients  have  been  proposed  (see  Cazacu  and  Barlat,  2003;  Barlat 
et  al.,  2005;  Hu,  2005;  etc).  For  metals  with  cubic  crystal  structure,  the  basic  deformation 
mechanism  is  slip.  Since  slip  does  not  depend  on  the  sign  of  the  shear  stress  i.e.  can  operate 
both  forward  and  backward,  the  tensile  and  compressive  yield  stresses  are  equal  so  the 
yield  functions  are  symmetric  about  the  origin  in  the  stress  space.  Metals  with  hexagonal 
close  packed  (hep)  crystal  structure  deform  plastically  by  slip  and  twinning.  As  opposed  to 
slip,  twinning  is  a  directional  shear  mechanism:  shear  in  one  direction  can  cause  twinning 
while  shear  in  the  opposite  direction  cannot.  For  example,  in  magnesium  alloys  sheets 
twinning  is  not  active  in  tension  along  any  direction  in  the  plane  of  the  sheet,  but  is  easily 
activated  in  compression.  As  a  result  the  average  initial  in-plane  compressive  yield  stress  is 
about  half  the  average  in-plane  tensile  yield  stress  (e.g.  see  Lou  et  al.,  2007).  Thus,  the  yield 
surfaces  are  not  symmetric  with  respect  to  the  stress  free  condition.  Since  hep  metals  sheets 
exhibit  strong  basal  textures  (c-axis  oriented  predominantly  perpendicular  to  the  thickness 
direction),  a  pronounced  anisotropy  in  yielding  is  observed.  To  account  for  both  strength 
differential  (SD)  effects  and  the  anisotropy  displayed  by  hep  metals,  Hosford  (1966)  pro¬ 
posed  the  following  modification  of  Hill’s  (1948)  orthotropic  yield  criterion: 

At 7 xx  +  BOyy  +  (-B  -  A)azz  +  F(<jyy  —  ffz-)2  +  G(azz  -  oxxf  +  H (axx  -<rn,)2  =  l,  ( 1 ) 

where  A,  B,  F,  G,  H  are  material  coefficients  and  x,  y,  z  are  normal  to  the  mutually  orthog¬ 
onal  planes  of  symmetry  of  the  material  .  Since  the  criterion  does  not  involve  shear  stres¬ 
ses,  it  cannot  account  for  the  continuous  variation  of  the  plastic  properties  between  the 
material  axes  of  symmetry.  Liu  et  al.  (1997)  have  proposed  an  extension  of  Hill  (1948) 
yield  criterion  in  the  form: 

{F(oyy  -  azz)2  +  G{azz  -  axxf  +  H (axx  -  u^)2  +  2Lo^  +  2 Ma2xz  +  2No]y }1/2 

+  lGXX  +  Jdyy  A  K(JZZ  =  1  .  (2) 

In  the  above  equation,  F,  G,  FI ,  L,  M,  N,  /,  /,  K  are  independent  material  coefficients. 
Although  this  yield  criterion  captures  SD  effects,  the  asymmetry  in  yielding  is  due  to  pres¬ 
sure  effects.  The  effects  of  hydrostatic  pressure  on  macroscopic  plastic  flow  have  been 
reported  by  Spitzig  and  Richmond  (1984)  for  fully  dense  metals  (e.g.  steels).  However, 
these  effects  are  small  at  low  pressure  levels  and  lead  to  yield  stresses  which  are  larger 
in  compression  than  in  tension.  For  HCP  metals,  the  predominant  mechanism  responsible 
for  SD  effects  is  twinning  (Hosford  and  Allen,  1973),  which  typically  leads  to  lower  initial 
yield  stresses  in  compression  than  in  tension  for  in-plane  loadings  of  rolled  sheets.  There¬ 
fore,  the  criterion  in  Eq.  (2)  may  not  capture  with  accuracy  the  behavior  of  HCP  metals. 
Recently,  based  on  results  of  polycrystalline  simulations,  Cazacu  and  Barlat  (2004)  have 
proposed  a  macroscopic  isotropic  yield  criterion  expressed  in  terms  of  the  invariants  of 
the  stress  deviator,  which  captures  the  asymmetry  in  yielding  between  tension  and  com¬ 
pression.  To  describe  both  the  asymmetry  and  anisotropy  in  yielding  of  magnesium  and 
magnesium  alloys  sheets,  an  extension  of  this  criterion  to  orthotropy  was  formulated  using 
the  generalized  invariants  approach  proposed  by  Cazacu  and  Barlat  (2001).  For  full  stress 
state  (3D)  conditions,  anisotropy  is  described  by  eighteen  coefficients.  This  orthotropic 
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yield  function  is  homogeneous  of  degree  three  in  stresses,  yet  for  certain  hep  materials  such 
as  titanium  the  yield  surfaces  are  quadratic  (see  Liu  et  ah,  1997;  Cazacu  et  ah,  2006;  etc.). 
To  overcome  this  limitation,  Cazacu  et  ah  (2006)  proposed  an  isotropic  yield  criterion  for 
which  the  degree  of  homogeneity  is  not  fixed.  To  capture  simultaneously  anisotropy  and 
tension/compression  asymmetry,  this  isotropic  yield  criterion  was  extended  to  orthotropy 
by  applying  a  fourth-order  linear  operator  on  the  stress  deviator.  Thus,  for  full  3-D  stress 
states,  nine  anisotropy  coefficients  are  involved  in  the  criterion. 

This  paper  presents  full  stress  3-D  yield  criteria  for  describing  the  anisotropic  plastic 
response  of  textured  metals.  The  aim  is  to  develop  models  that  can  be  applicable  to  mate¬ 
rials  that  exhibit  strength  differential  effects  as  well  as  to  materials  for  which  not  noticeable 
difference  exist  between  the  behavior  in  tension  and  compression  under  monotonic  load¬ 
ing.  Key  in  this  development  is  the  use  of  the  isotropic  yield  function  of  Cazacu  et  ah 
(2006).  Anisotropy  is  introduced  using  several  linear  transformations.  In  the  next  section, 
this  isotropic  yield  criterion  is  succinctly  presented.  After  reviewing  general  aspects  of  lin¬ 
ear  transformations  operating  on  the  Cauchy  stress  tensor,  in  Section  3,  a  new  anisotropic 
yield  function  involving  two  linear  transformations  is  introduced.  The  input  data  needed 
for  the  calculation  of  anisotropic  yield  function  coefficients  are  discussed.  Illustrative 
examples  of  application  of  this  yield  function  to  the  description  of  anisotropy  and  ten¬ 
sion/compression  asymmetry  of  materials  with  cubic  and  hexagonal  close  packed  crystal 
structure  are  presented.  It  is  shown  that  by  incorporating  more  than  two  linear  transfor¬ 
mations  in  the  isotropic  yield  criterion,  a  very  high  degree  of  accuracy  in  the  representa¬ 
tion  of  a  very  large  set  of  anisotropic  tensile  and  compressive  data  can  be  achieved. 
Moreover,  it  is  shown  that  the  3-D  yield  criterion  involving  two  linear  transformations 
captures  with  high  accuracy  the  anisotropy  in  yielding  of  metals  with  no  tension/compres¬ 
sion  asymmetry. 

2.  Isotropic  yield  criterion 

According  to  classical  results  of  the  theory  of  representation  (e.g.  Wang,  1970),  the 
yield  function  of  an  isotropic  pressure-insensitive  material  must  be  expressed  in  terms  of 
the  three  invariants  of  S,  the  deviator  of  stress  tensor  a.  Many  sets  of  three  invariants 
are  possible  (Zyczkowski,  1981)  but,  in  this  work,  the  principal  values  Sk  are  used  to  rep¬ 
resent  the  isotropic  yield  function  /  i.e. 

f(S)=F(SuS2,Si),  (3) 

where  Sh  i  =  \  ...  3  are  the  principal  values  of  the  stress  deviator  S.  I-  must  be  an  isotropic 
function  of  its  arguments,  i.e., 

F(Si,Sn,Sm)  =  F(Si,S2,Sj), 

where  (1, 11,111)  are  permutations  of  (1,2,3). 

If  a  material  only  deforms  by  a  reversible  shear  mechanism  such  as  slip,  yielding 
depends  only  on  the  magnitude  of  the  resolved  shear  stress.  Thus,  f(S)  =  /(  —S)  and  the 
yield  locus  in  the  deviatoric  n  plane  (plane  which  passes  through  the  origin  and  is  perpen¬ 
dicular  to  the  hydrostatic  axis)  must  have  sixfold  symmetry.  If  the  material  deforms  by 
twinning  yielding  depends  on  the  sign  of  the  applied  shear  stress,  i.e.,  yield  in  tension 
and  compression  should  be  different.  Hosford  and  Allen  (1973)  used  a  modified  Taylor 
polycrystal  model  to  show  that  for  randomly  oriented  fee  and  bcc  crystals  deforming 
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solely  by  twinning,  the  ratio  between  the  yield  stress  in  tension  and  compression  should  be 
0.78  and  1.28,  respectively.  To  account  for  this  strength  differential  effect,  the  following 
isotropic  yield  function  was  proposed  by  Cazacu  et  al.  (2006) 

0  =  (N  -  kS\)a  +  (|S2|  -  kS2)a  +  (|S3|  -  kS3)a,  (4) 

where  a  and  k  are  material  parameters.  It  was  shown  that  for  a  fixed  value  of  the  degree  of 
homogeneity  a ,  the  parameter  k  is  expressible  solely  in  terms  of  the  ratio  between  erT,  the 
uniaxial  yield  in  tension,  and  <rc  the  uniaxial  yield  in  compression,  respectively,  i.e. 


/  2a-2(aT/ac)a\ ° 
\  (2-aT / ocY -2  j 

f  2a -2-(ax  /  QcY  1  ^ 
(  (2  ot/ocY-2  J 


(5) 


It  was  shown  that  this  yield  function  is  also  in  excellent  agreement  with  results  obtained 
using  the  self-consistent  polycrystal  (VPSC)  model  of  Lebensohn  and  Tome  (see  for  exam¬ 
ple  Lebensohn  and  Tome,  1993)  for  randomly  oriented  polycrystals  deforming  solely  by 
twinning.  Assuming  a  =  2,  and  using  (5),  the  strength  differential  parameter  is 
k  =  —  0.31  for  fee  materials,  while  k  =  0.3 1  for  bcc  materials.  Furthermore,  for  a  =  3 
and  k  —  —0.0645,  the  yield  loci  (4)  are  in  excellent  agreement  with  the  yield  loci  for  ran¬ 
domly  oriented  Zr  (hep)  polycrystals  deforming  solely  by  tensile  twinning 
{ 1 0 1 2}  ( 1 0 1  1 )  and  compressive  twinning  {112  2}  (112  3)  obtained  using  the  VPSC 
model. 

Irrespective  of  the  value  of  a,  if  the  yield  stresses  in  tension  and  compression  are  equal 
k  =  0.  In  particular,  for  k  =  0  and  a  =  2,  the  yield  criterion  (4)  reduces  to  the  Von  Mises 
yield  criterion. 


3.  Anisotropic  yield  criteria 


To  extend  this  isotropic  criterion  to  orthotropy,  Cazacu  et  al.  (2006)  applied  a  linear 
transformation  on  the  deviatoric  stress  tensor  S,  i.e.  in  Eq.  (4)  Si,  S2,  S3  were  substituted 
by  the  principal  values  of  a  transformed  tensor  £  defined  as 

£  =  C:S.  (6) 

Thus,  the  orthotropic  yield  criterion  (denoted  in  the  following  as  CPB06)  is  of  the  form: 

F  =  <i>(zuz2,z3)  =  (1^1  -  kzxy  +  (|z2|  -  kz2y  +  (|S3|  -  ki3y,  (7) 

where  Xu  Z2,  S3  are  the  principal  values  of  X.  The  only  restrictions  imposed  on  the  fourth- 
order  tensor  C  (see  Eq.  (6))  are:  (i)  to  satisfy  the  major  and  minor  symmetries  and  (ii)  to  be 
invariant  with  respect  to  the  orthotropy  group.  Thus,  for  3-D  stress  conditions  CPB06  in¬ 
volves  nine  independent  anisotropy  coefficients  and  it  reduces  to  the  isotropic  criterion  (4) 
when  C  is  equal  to  the  fourth-order  identity  tensor.  It  is  worth  noting  that  although  the 
transformed  tensor  is  not  deviatoric,  the  orthotropic  criterion  is  insensitive  to  hydrostatic 
pressure  and  thus  the  condition  of  plastic  incompressibility  is  satisfied  (see  Cazacu  et  ah, 
2006  for  details).  For  k  G.  [—1,1]  and  any  integer  a  >  1,  the  anisotropic  yield  function  is 
convex  in  the  variables  Z\,  X2,  T3.  The  experimental  data  on  hep  materials  reported  in 
the  literature  (e.g.  Lee  and  Backofen,  1966;  Kelley  and  Hosford,  1968)  is  generally  limited 
to  compressive  and  tensile  yield  stresses  along  the  axes  of  symmetry  of  the  respective  mate¬ 
rials.  It  was  shown  (see  Cazacu  et  al.,  2006)  that  the  CPB06  yield  criterion  describes  with 
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accuracy  the  tension/compression  asymmetry  of  these  materials  as  well  as  that  of  high- 
purity  zirconium  with  in-plane  isotropy  (Plunkett  et  al.,  2006).  However,  it  is  clear  that 
nine  anisotropy  coefficients  may  not  be  sufficient  to  describe  with  very  good  accuracy 
yielding  of  materials  with  very  pronounced  in-plane  anisotropy. 

To  increase  the  number  of  anisotropy  coefficients  in  the  formulation,  instead  of  one  linear 
transformation,  n  linear  transformations  (n  ^  2)  can  be  considered.  Such  a  methodology 
was  used  by  Barlat  et  al.  (2003,  2005,  2006)  and  Bron  and  Besson  (2004)  to  describe  the  pro¬ 
nounced  anisotropy  displayed  by  certain  aluminum  alloys.  For  example,  Yld2004-18p 
proposed  by  Barlat  et  al.  (2005)  involves  two  linear  transformations.  When  the  two  transfor¬ 
mations  are  equal,  Yld2004-18p  reduces  to  Yld91  orthotropic  yield  criterion  (Barlat  et  al., 
1991).  Here,  we  demonstrate  that  additional  linear  transformations  can  be  incorporated  into 
the  CPB06  criterion  for  an  improved  representation  of  the  anisotropy  yield  surface. 

The  following  analytic  yield  function,  denoted  CPB06ex2,  is  proposed: 

F(L, s')  =  (isj I  -  kXi)a  +  (|t2|  -  kz2y  +  (|s3|  -  kz}y  +  (\z\  \  -  k'z\y 

+  (\z'2\-kiz'2y  +  (\zi3\-k'z'3y,  (8) 


where  k  and  k'  are  material  parameters  that  allow  for  the  description  of  strength  differen¬ 
tial  effects,  a  is  the  degree  of  homogeneity,  while 

£  =  C:S  and  S'  =  C'  :  S.  (9) 


Let  ( x,y,z )  be  the  reference  frame  associated  with  orthotropy.  In  the  case  of  a  sheet,  jc, 
y,  and  z  represent  the  rolling,  transverse,  and  the  normal  directions.  Relative  to  the  ortho¬ 
tropy  axes  (x,y,  z),  the  fourth-order  tensors  C  and  C'  operating  on  the  stress  deviator  are 
represented  by 
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c' 
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0 

c' 

*—66 

(10) 


Thus,  for  3-D  stress  conditions  the  orthotropic  criterion  (8)  involves  18  anisotropy  coef¬ 
ficients.  When  Cu  =  1  and  all  other  CVj  =  0  (i  /  j)  and  k  =  k' ,  this  yield  function  reduces 
to  the  isotropic  yield  function  (4).  Note  that  when  C  =  C1  and  k  =  k'  the  proposed  crite¬ 
rion  reduces  to  the  CPB06  yield  criterion  (see  Eq.  (7)). 

Consider  uniaxial  loading  in  the  plane  (jc,j)  of  the  sheet  along  a  direction  at  angle  9 
with  the  rolling  direction  and  denote  by  erj  and  (7g  the  tensile  and  compressive  yield  stres¬ 
ses,  respectively.  Then,  according  to  the  proposed  criterion  (8): 
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[Mil  -  +  [M2|  -  kA2]a  +  [\A3\  -  kA3]a 

+[Mil  -  k'A\y  +  [Mil  -  k'A'2y  +  [Mil  -  £'4r 


w  I  [\Ax\  +  kAx}a  +  [\A2\  +  kA2]a  +  [Ms |  +  kA3]a  r 

(  +[Mi I  +  +  [Mil  +  ^i]°  +  [Mil  +  k'AY ) 

where  XT  is  the  tensile  yield  stress  in  the  rolling  direction  (i.e.  for  6  =  0),  and 


(11) 


(12) 


At  —  ^  Mh  +  Ag  +  \J (Axx  —  Ayy)  +  4Axl, J  ,  A2  —  2  —  \J (Axx  —  Aly)  +  4AX^J ,  A3  —  A„, 

Axx  =  <P  1  cos2  0+  'Pi  sin2  0,  =  <P2  cos2  0  +  sin2  0,  ^4z;  =  <f>3  cos2  0  +  P3  sin2  0, 

=  Cee  sin  0  cos  0,  and 

2  11  2  11  2  11 

j  C32  —  j  Cu  —  ^€13,  P2  =  j  C22  —  j  C12  —  j  C23,  ^3  =  j  C23  —  j  C13  — -C33. 

(13) 


Relations  similar  to  (13)  express  <P\ .  </>(,.  in  terms  of  the  anisotropy  coefficients  C\j  and 
A\,A'2,A'3  as  a  function  of  the  angle  0  and  the  anisotropy  coefficients  C1-,  respectively. 

In  particular,  X°,  the  compressive  yield  stress  in  the  rolling  direction  is 

1 

[|<Z>i|+  k<l>{\a  +  [|<P2|  +  k<P2]a  +  +  k<P3]a 

+  [KI  +  k'$[]a  +  [\&2\  +  k'<P'2]a  +  [|<2>'|  +  k'$'3]a 


where 


A>\  —  —  Cu  —  C 1 2  —  -  C|3,  —  —  d7 12  —  37  C22  —  t-  C23,  @3  —  —  C2 13  —  —  C23  —  —  C733 . 

3  3  3’  -  3  3  3’  3  3  3 

(15) 

If  Yt  and  yx  are  the  tensile  and  compressive  yield  stresses  in  the  j-direction  (transverse 
direction)  then 


yT  |  [i-Pii-^ir+n^i-^r+o^i-^r 

{  +[\r1\-k'rl}a  +  [\r2\-k'r2}a  +  [\'p'3\-k,r3\a 

yC  =  xT ]  ph|T WiY  +  [| t2\  +  +  [|  <f3|  +  kv3 f 

{  +[!¥",!  +  k!rx\a  +  [\r2\  +  k'r2]a  +  [|  r3\  +  k'r3)a 

If  ZT  and  ZF  are  the  tensile  and  compressive  yield  stresses  in  the  "-direction  (normal  to  the 
sheet)  then 
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j _ I _ 

]  [|77|  |  -  A-77,r  +  [|772|  -  hTI2]a  +  [|i73|  -  kll3)a 

{  +[\n\\- k'n[]a  +  [\n'2\  - k'n'2]a  +  [\n'3\  - k'n'3]a 


i  [|^i|+^^i]  ~k^n2\+kn2\  +  [|773|  +  AT73] 

l  +[\n\  I  +  k'n\r  +  [|27'|  +  k'n'2]°  +  [\n'3\  +  k'n’Y 


where 


III  —  (*1 3  C]  |  —  —  Cp, 
3  3  3 


III  —  Ci3  —  —  C12  ~  —  Ci 2,  III  —  :r  C33  —  —  C13 

3  3  3  3  3 


(18) 


(19) 


1 

3 


C23 


(20) 


and  relations  similar  to  (20)  express  each  77'  term  as  a  function  the  anisotropy  coefficients 
C\j.  Due  to  plastic  incompressibility,  yielding  under  equibiaxial  tension  and  compression 
occurs  when  axx  and  avy  are  both  equal  to  Zc  and  ZT,  respectively. 

The  yield  stress  in  pure  shear  in  the  sheet  plane  is 

{ [I ^66 1  +  kC66}“  +  [|C66|  -  kC66]a  +  [|Cg6|  +  kC'66]a  +  [|C'6|  -  k'C'm]°} 

(21) 

and  the  yield  stresses  in  shear  with  respect  to  the  other  directions  of  orthotropy  can  be 
found  by  the  substitution  of  C44,  and  C'44  or  C55,  and  C'55  in  the  place  of  C66,  and  C'66. 

Furthermore,  we  assume  that  the  plastic  potential  coincides  with  the  yield  function.  Let 
denote  by  r0  the  Lankford  coefficient  (width  to  thickness  strain  ratios)  under  uniaxial  ten¬ 
sile  or  compressive  loading  in  a  direction  at  angle  0  with  the  rolling  direction  in  the  plane 
(xy).  Hence 


sin2  sin(20)  +  cos2  9 

OGxX  V  7  OGXy  OGyy 

Si7  |  dF' 

dG XX  '  dGyy 


(22) 


where  is  given  by  Eq.  (8)  and  the  explicit  expressions  of  the  derivatives  in  terms  of 

the  components  of  the  Cauchy  stress  are  given  in  Appendix  A. 

Similarly,  more  linear  transformations  can  be  incorporated  in  the  orthotropic  yield  cri¬ 
terion  (8)  (denoted  CPB06ex3  for  3  transformations,  etc.)  until  a  desired  level  of  accuracy 
in  the  description  of  the  plastic  anisotropy  is  obtained,  i.e. 

T’(S)  =/(1)+/(2)  +  ---+/(p)+/W,  withl^p^n,  (23) 

where /w  =  (|T^|  -  k^Z^Y  +  (|z'rt|  -  k^zff  +  (|zf  |  -  k^zff.  In  Eq.  (23)  k(p)  are 
material  parameters,  z[p) ,i=  1  ...  3  are  the  principal  values  of  =  C^’S  with  C(p)  being 
the  fourth-order  orthotropic  tensor  associated  to  the  pih  linear  transformation. 

The  anisotropy  coefficients  involved  in  the  yield  criteria  can  be  found  through  the  min¬ 
imization  of  an  error  function.  The  experimental  data  in  the  error  function  may  consist  of 
flow  stresses  and  r- values  in  tension  and  compression  corresponding  to  different  orienta¬ 
tions  in  the  plane  of  the  sheet,  biaxial  flow  stress  in  tension  and  compression,  as  well  as 
out-of-plane  yield  stresses  and  /--values.  If  experimental  data  are  not  available  for  a  given 
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strain  path,  they  can  be  substituted  with  numerical  data  obtained  based  on  polycrystalline 
calculations 


2 


(24) 


In  the  above  equation  i  represents  the  number  of  experimental  yield  stresses,  j  represents 
the  number  of  experimental  /'-values  available  while  the  superscript  indicates  whether  the 
corresponding  value  is  experimental  or  predicted. 

In  the  next  section,  we  will  apply  the  CPB06ex2  CPB06ex3,  and  CPB06ex4  formula¬ 
tions  to  the  description  of  the  anisotropy  and  tension  compression  asymmetry  of  HCP, 
BCC,  and  FCC  alloys. 

4.  Application  to  materials  exhibiting  asymmetry  between  tension  and  compression 

4.1.  AZ3 IB  magnesium  (HCP) 

Lou  et  al.  (2007),  reported  results  of  an  experimental  investigation  of  the  monotonic 
and  cyclic  mechanical  behavior  of  a  commercial  AZ3 1 B  magnesium  alloy  sheet  in  tension, 
compression,  and  simple  shear.  The  characterization  of  the  in-plane  anisotropy  of  the 
sheet  was  done  by  performing  uniaxial  tensile  and  compression  tests  in  the  rolling  direc¬ 
tion,  at  45°  to  the  rolling  direction  and  in  the  transverse  direction.  The  material  displays 
significant  anisotropy  as  well  as  asymmetric  yield  and  hardening  behavior  which  is  asso¬ 
ciated  with  the  polarity  of  twinning  (the  easily  activated  twin  { 1 0 1 2}  is  a  tensile  twin). 
In  this  paper,  we  are  only  concerned  with  modeling  with  accuracy  the  initial  yielding  of 
the  material.  The  data  reported  at  0.2%  olfset  are  sufficient  for  constructing  the  ten¬ 
sion-tension  and  compression-compression  quadrants  of  the  experimental  plane-stress 
yield  locus.  We  will  now  compare  Lou  et  al.  (2007)  to  the  theoretical  yield  locus  calculated 
with  the  proposed  orthotropic  yield  criterion  CPB06ex2  (Eq.  (8))  in  Fig.  1.  The  parameters 
involved  in  the  expression  of  the  theoretical  yield  loci  for  in-plane  (xy)  stress  states  were 
determined  using  equations  (1.11)  and  (1.12)  and  all  the  available  data.  For  this  material, 
a  =  4.  The  values  of  these  anisotropy  coefficients  are  given  in  Table  1.  Fig.  2  compares 
experimental  and  calculated  yield  stresses  normalized  by  the  uniaxial  tensile  stress  in  the 
rolling  direction  (0.2%  offset)  and  r- values  as  a  function  of  the  loading  direction  (angle 
from  the  RD).  Since  all  available  experimental  data  were  used  to  identify  the  coefficients 
involved  in  the  criterion,  the  calculations  are  not  really  predictive.  However,  it  is  clearly 
seen  that  the  CPB06ex2  yield  criterion  is  flexible  enough  to  describe  the  whole  set  of  aniso¬ 
tropic  data  (12  independent  material  properties,  i.e.  both  /--values  and  yield  stresses)  for 
both  tensile  and  compressive  loadings.  This  flexibility  cannot  be  achieved  with  models  that 
are  based  on  one  linear  transformation. 

4.2.  Medium  carbon  low  alloy  steel  (BCC) 

The  second  application  example  of  CPB06ex2  yield  criterion  (see  Eq.  (8))  pertains  to  a 
medium  carbon  alloy  steel  plate  for  which  Benzerga  et  al.  (2004)  reported  asymmetry 
between  the  yield  behavior  in  tension  and  compression.  The  anisotropy  of  deformation 
was  characterized  through  a  series  of  tension  and  compression  tests  along  the  rolling 
(x),  transverse  (>’).  and  normal  direction  (z)  as  well  as  along  45°  orientations  in  the  (xy), 
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Fig.  1.  Plane-stress  yield  loci  (axy  =  0)  for  AZ31B  Mg  alloy  according  with  the  proposed  CPB06ex2  yield 
criterion  and  experiments  (symbols).  Data  after  Lou  et  al.  (2007). 

(yz),  and  (zx)  planes,  respectively.  Under  off-axes  loading  in  any  symmetry  plane,  the 
strain  ratios  were  defined  as  the  in-plane  radial  strain  over  the  out-of-plane  radial  strain. 
For  this  BCC  crystal  structure  material,  it  was  observed  that  in  compression  there  is  no 
significant  difference  between  the  yield  stresses  along  the  x,  y,  and  z  orientations  although 
the  corresponding  /--values  are  different  (see  Tables  3  and  4).  All  the  available  experimental 
data  ( 14  measured  yield  stress  and  /'-values)  were  used  to  determine  the  coefficients  of  the 
CP05ex2  yield  criterion.  For  this  material,  a  =  5.  The  values  of  these  anisotropy  coeffi¬ 
cients  are  given  in  Table  2. 

Fig.  3  shows  the  theoretical  yield  surface  in  the  normalized  stress  plane  (ax,  ay)  along 
with  the  experimental  data  points.  In  Tables  3  and  4  are  given  the  values  of  the  calculated 
and  measured  yield  stresses  and  Lankford  coefficients.  Note  that  the  proposed  theory  cap¬ 
tures  very  well  the  observed  asymmetry  in  yielding  between  tension  and  compression,  the 
strong  anisotropy  in  /--values  and  a  relatively  mild  anisotropy  in  yield  stresses  for  both  ten¬ 
sion  and  compression.  This  flexibility  cannot  be  achieved  with  Hill’s  (1948)  yield  criterion 
which  is  generally  used  for  modeling  steels. 

4.3.  2090-T3  aluminum  (FCC) 

It  was  reported  by  Barlat  et  al.  (1991)  that  for  a  cold  rolled  2090-T3  aluminum  alloy 
sheet  under  biaxial  plane-stress  loading  conditions,  the  tension-tension  (I)  and  compres- 


Table  1 

CPB06ex2  yield  criterion  parameters  for  AZ31B  Mg  alloy 


k 

Cn 

C\2 

^*13 

C22 

Q23 

C33 

c66 

0.295 

0.771 

— 0.914 

— 1.312 

0.654 

-1.209 

-0.293 

1.617 

kf 

c' 

C11 

c' 

'“'12 

C' 

'-'13 

c' 

'“'22 

c' 

'"'23 

c' 

*-33 

c' 

u66 

0.900 

0.771 

0.337 

— 0.764 

0.646 

— 0.401 

-0.022 

0.499 
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Fig.  2.  Anisotropy  of  the  yield  stress  (normalized  by  the  tensile  stress  in  the  rolling  direction]  and  the  r-values  for 
AZ31B  Mg,  measured  and  predicted  with  the  CPB06ex2  yield  criterion  for  (a)  tensile  loading;  (b)  compressive 
loading.  Experimental  data  after  Lou  et  al.  (2007). 


sion-compression  quadrants  (IV)  are  not  identical  due  to  a  yielding  difference  in  tension 
and  compression.  However,  this  SD  effect  is  a  snapshot  of  the  material  behavior  at  yield 
due  to  micro-plasticity  effects  related  to  prior  thermo-mechanical  processing.  Aluminum 
alloys  deformed  beyond  1-2%  are  expected  to  behave  the  same  (identical  stress-strain 
curves)  in  tension  and  compression.  The  deformation  mechanisms  in  tension  and  compres¬ 
sion  are  not  fundamentally  different  in  aluminum  alloys,  in  contrast  with  HCP  materials. 
This  example  is  provided  in  order  to  show  the  capability  of  the  model  to  capture  simulta¬ 
neously  the  SD  as  well  as  strong  anisotropy  effects  for  a  material  that  was  characterized 
extensively,  i.e.,  with  a  large  number  of  experimental  data  to  account  for. 


Table  2 

CPB06ex2  coefficients  for  a  medium  carbon  low  alloy  steel 


k 

c„ 

Cl2 

Cl3 

C22 

^23 

C33 

C44 

C55 

0,6 

0.030 

1.283 

0.208 

-0.421 

0.535 

-0.791 

0.746 

1.418 

1.329 

1.156 

V 

r' 

C11 

c' 

12 

c' 

u13 

c' 

'“'22 

c' 

'“'23 

c' 

C33 

c' 

*“44 

c' 

u55 

r< 

c66 

-0.027 

1.283 

0.380 

-0.001 

2.145 

0.792 

0.218 

1.418 

1.329 

1.721 
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Fig.  3.  Projection  of  the  CPB06ex2  yield  surface  in  the  (<rv,oy)  plane  for  a  medium  carbon  low  alloy  steel. 
Symbols  represent  normalized  experimental  data  of  Benzerga  et  al.  (2004). 


Table  3 

Measured  and  theoretical  normalized  yield  stress  values  for  a  medium  carbon  low  alloy  steel.  Experimental  data 
after  Benzerga  et  al.  (2004) 


Tension 

Compression 

XT 

yt 

A'c 

>c 

Experiment 

1 

0.914 

0.963 

0.960 

0.958 

CPB06ex2 

1 

0.914 

0.963 

0.960 

0.958 

Table  4 

Measured  and  predicted  /--values  for  a  medium  carbon  low  alloy  steel 


Loading  direction 

Tension 

Compression 

x  y 

y  z 

45°in  (xy)  plane 

45°in  (x~)  plane 

45°in  (yz)  plane 

Experiment 

0.67  1.4 

0.71 

1.3  1.05 

0.91 

1.45 

1.3 

CPB06ex2 

0.67  1.4 

0.71 

1.3  1.05 

0.91 

1.45 

1.3 

Experimental  data  after  Benzerga  et  al.  (2004). 


The  experimental  data  consists  of  tensile  and  compressive  yield  stresses  measured  along 
seven  directions  in  the  plane  of  the  sheet  (i.e.  from  rolling  to  transverse  directions  in  15° 
increments).  The  strain-ratio  anisotropy  was  measured  only  for  tensile  conditions.  Biaxial 
flow  stress  from  a  bulge  test,  and  the  biaxial  r  value  from  the  disk  compression  test  (rh) 
were  reported  in  Barlat  et  al.  (2003).  The  description  of  the  uniaxial  tension  properties 
of  this  material  to  a  desired  degree  of  accuracy  has  posed  challenges.  Many  anisotropic 
yield  functions  of  different  levels  of  complexity  have  been  proposed  in  order  to  capture 
both  the  anisotropy  in  the  tensile  flow  stresses  and  tensile  /'-values.  However,  with  the 
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exception  of  the  yield  function  proposed  by  Bron  and  Besson  (2004)  and  Barlat  et  al. 
(2005)  most  of  the  formulations  were  developed  only  for  plane-stress  states  or  do  not  fulfill 
automatically  the  requirements  for  convexity.  To  describe  the  anisotropic  tensile  proper¬ 
ties  of  the  material  and  at  the  same  time  account  for  the  difference  in  flow  stresses  between 
tension  and  compression  in  the  rolling  and  transverse  directions,  Yoon  et  al.  (2000)  used 
the  anisotropic  yield  function  Yld96  (Barlat  et  al.,  1997)  and  kinematic  hardening.  How¬ 
ever,  with  this  approach,  the  accuracy  in  representing  the  tensile  and  compressive  flow 
stresses  for  off-axes  loading  in  the  plane  of  the  sheet  (i.e.  in  directions  other  than  the 
RD  and  TD)  is  reduced.  In  this  paper,  we  will  show  that  by  introducing  several  linear 
transformations  in  the  isotropic  yield  criterion  (4)  (denoted  CPB06ex3  for  three  linear 
transformations,  CPB06ex4  for  four  linear  transformations)  the  desired  level  of  accuracy 
in  the  description  of  the  plastic  anisotropy  for  both  tensile  flow  stresses,  compressive  flow 
stresses  and  tensile  r- values  can  be  achieved.  Specifically,  we  model  the  material  using  the 
proposed  yield  criterion  with  2,  3  or  4  transformations,  respectively.  For  the  determination 
of  the  material  parameters  involved  in  each  yield  criterion  we  use  all  the  available  exper¬ 
imental  data  (a  total  of  22  data  points).  For  this  material,  the  exponent  a  =  12.  A  high 
value  for  the  homogeneity  exponent  was  necessary  to  capture  with  accuracy  the  material’s 
response.  The  value  of  the  exponent  “a”  was  adjusted  until  the  best  approximation  of  the 
entire  data  set  was  obtained.  For  each  yield  criterion  CPB06ex2,  CPB06ex3,  CPB06ex4 
and  the  values  of  the  coefficients  associated  with  strength  differential  effects  are  listed  in 
Tables  5-7,  respectively. 

Fig.  4  shows  the  plane-stress  yield  loci  (xxy  —  0)  calculated  using  CPB06ex2  (14  anisot¬ 
ropy  coefficients  for  plane  stress),  CPB06ex3  (21  anisotropy  coefficients  for  plane  stress) 
and  CPB06ex4  (28  anisotropy  coefficients  for  plane-stress  states)  along  with  the  experi¬ 
mental  flow  stresses.  Fig.  5  shows  the  anisotropy  of  the  tensile  and  compressive  flow  stres¬ 
ses  and  tensile  r-values  as  described  by  CPB06ex2,  CPB06ex3,  and  CPB06ex4  in 
comparison  with  the  experimental  values.  It  can  be  noted  that  although  the  tension-com¬ 
pression  anisotropy  of  the  yield  loci  is  captured  with  accuracy  by  all  the  anisotropic  yield 
criteria  (see  Fig.  4),  it  is  clearly  seen  that  as  more  linear  transformations  are  involved  in  the 


Table  5 

CPB06ex2  coefficients  for  2090-T3 


k 

C,  i 

C\2 

Cl3 

C22 

^23 

C33 

0,6 

0.054 

0.453 

-0.841 

-1.248 

-1.058 

-2.284 

-3.201 

1.026 

K 

c' 

Hi 

c' 

'-'12 

c' 

c13 

c' 

'“'22 

c' 

'“'23 

c' 

'-'33 

c' 

'“'66 

0.027 

0.453 

-0.705 

1.148 

0.139 

-0.519 

0.878 

1.978 

Table  6 

CPB06ex3  coefficients  for  2090-T3 


k 

C„ 

C\2 

C13 

C22 

C23 

C33 

06 

-0.076 

4.192 

5.311 

4.031 

4.488 

5.443 

5.671 

2.074 

k! 

c 

Hi 

r' 

H2 

c' 

13 

c' 

'“'22 

c' 

'“'23 

c' 

'“'33 

c 

c66 

-0.020 

4.192 

5.095 

6.135 

6.750 

5.421 

5.696 

1.000 

k" 

c" 

Hi 

C" 

H2 

c" 

'-'  13 

c" 

'“'22 

r" 

'“'23 

c" 

u33 

c" 

c66 

-0.029 

4.192 

3.412 

2.562 

4.629 

2.794 

3.205 

1.422 
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Table  7 

CPB06ex4  coefficients  for  2090-T3 


k 

Cii 

C\2 

Cn 

C22 

C23 

Q3 

c66 

0.331 

0.417 

-0.709 

-0.588 

0.277 

-1.053 

-1.949 

1.275 

IP 

c' 

U11 

c' 

12 

^13 

c' 

'-'22 

c' 

'-'23 

c' 

*-'33 

c' 

c66 

-0.311 

0.417 

-0.442 

-0.549 

-0.975 

-1.793 

-0.924 

-0.820 

k" 

r" 

C" 

c" 

c13 

r" 

'-'22 

r" 

^23 

r" 

u33 

r" 

c66 

0.063 

0.417 

0.214 

1.295 

-0.057 

0.065 

-0.624 

1.924 

IP" 

ft III 

'-'ll 

ftlll 

*“  12 

/-till 

*■'13 

f till 
'-'22 

/—till 

'-'23 

ft  III 

*-'33 

ftlll 

c66 

0.625 

0.417 

-0.297 

0.856 

-0.570 

-0.160 

0.775 

0.938 

Fig.  4.  Comparison  between  experimental  and  theoretical  plane-stress  yield  surfaced  according  to  the  yield 
criteria  CPB06ex2,  CPB06ex3,  and  CPB06ex4  for  2090-T3  aluminum  alloy  sheet  (data  after  Yoon  et  al.,  2000). 


formulation  (i.e.  as  the  number  of  anisotropy  coefficients  is  increased)  the  accuracy  is  sig¬ 
nificantly  improved.  When  four  linear  transformations  are  introduced,  the  error  associated 
with  representing  all  22  experimental  data  points  tends  towards  zero. 

5.  Application  of  the  CPB06ex2  yield  criterion  to  textured  metal  sheets  exhibiting  symmetry 
between  tensile  and  compressive  properties 

Although  the  CPB06  yield  criterion  (Cazacu  et  al.,  2006)  and  the  extensions  proposed  in 
this  paper  were  formulated  to  capture  the  strength  differential  effects  most  often  associated 
with  materials  that  display  tension/compression  asymmetry  in  yielding,  the  proposed  for¬ 
mulations  are  not  limited  to  such  materials.  If  the  yield  in  tension  is  equal  to  the  yield  in 
compression,  the  parameters  k  and  k'  associated  with  strength  differential  effects  are  auto¬ 
matically  zero. 
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Fig.  5.  Anisotropy  of  both  tensile  and  compressive  flow  stresses  (normalized  by  uniaxial  tension  in  the  rolling 
direction!  and  tensile  /--values  for  2090-T3  alloy  sheet,  measured  (data  after  Yoon  et  al.,  2000)  and  predicted  by 
the  yield  criteria  CPB06ex2,  CPB06ex3  and  CPB06ex4. 


Ill  this  section,  we  apply  the  proposed  CPB06ex2  yield  criterion  (Eq.  (8))  to  materials 
with  cubic  crystal  structure  (FCC  and  BCC)  for  which  it  is  commonly  assumed  that  the 
tensile  and  compressive  yield  stresses  are  equal.  The  first  example  pertains  to  61 1 1-T4  alu- 


Table  8 

CPB06ex2  yield  criterion  parameters  for  6111-T4  Al  alloy 


Cn 

C\2 

Cl3 

Cl2 

C23 

C33 

C44 

^55 

c66 

0.250 

-1.399 

-0.971 

-0.375 

-2.254 

-0.831 

1.403 

1.101 

1.316 

c' 

C11 

c' 

u12 

c' 

u13 

c' 

^22 

c' 

'-'23 

c' 

u33 

c' 

^44 

c' 

c55 

c' 

c66 

0.250 

-0.840 

-1.330 

0.780 

-0.565 

0.479 

1.094 

1.448 

1.413 
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minum  alloy  sheet  sample  (data  after  Barlat  et  al.,  2005).  Tensile  yield  stresses  and  r-values 
along  seven  directions  in  the  plane  of  the  sheet  (i.e.  from  rolling  to  transverse  directions  in 
15°  increments)  were  measured  from  uniaxial  tensile  tests.  Additionally,  the  experimental 
values  of  the  balanced  biaxial  yield  stress  as  well  as  the  corresponding  /'-value  obtained 
using  a  disk  compression  test  were  reported.  In  addition  to  the  experimental  data,  VPSC 
polycrystal  calculations  were  carried  out  to  determine  uniaxial  yield  stresses  at  45°  and  the 
yield  stresses  in  simple  shear  in  the  Y-Z  and  X  Z  planes.  For  this  material  the  exponent 
a  =  12,  while  the  values  of  the  normalized  anisotropy  coefficients  are  given  in  Table  8. 
Table  9  lists  the  reported  VPSC  yield  stresses  and  corresponding  CPB06ex2  predictions. 
Figs.  6  and  7  compare  the  CPB06ex2  plane-stress  yield  locus  (axy  =  0)  and  the  predicted 
uniaxial  flow  properties  with  the  experimental  data,  respectively.  The  measured  rb-value 
is  1.225,  while  the  CPB06ex2  yield  criterion  prediction  is  of  1.227.  The  very  good  agree¬ 
ment  shows  that  CPB06ex2  is  flexible  enough  to  describe  the  whole  set  of  anisotropic  data 
(r-values  and  yield  stresses)  simultaneously. 

In  the  previous  section,  we  modeled  the  strongly  textured  2090-T3  aluminum  alloy  sam¬ 
ple  taking  into  consideration  the  observed  asymmetry  between  tension  and  compression  of 
this  material.  However,  in  the  literature,  this  tension/compression  asymmetry  is  generally 
neglected.  Next,  we  assume  that  the  yield  surface  of  the  material  is  symmetric  about  the 
origin  in  the  stress  space  and  apply  the  yield  criterion  CPB06ex2  with  k  =  0  and  k  =  0. 
The  values  of  the  coefficients  of  anisotropy  for  this  material  are  given  in  Table  10 


Table  9 

Comparison  between  theoretical  normalized  yield  stress  values  for  6111-T4  Al  alloy 


Y-Z  45°  tension 

Z-X  45°  tension 

Y-Z  simple  shear 

Z-X  simple  shear 

VPSC 

1.07 

1.04 

0.67 

0.65 

CPB06ex2 

1.07 

1.04 

0.67 

0.65 

VPSC  results  after  Barlat  et  al.  (2005). 


Fig.  6.  CPB06ex2  (k  =  0,17  =  0)  yield  surfaces  for  plane  stress  in  the  X  Y  plane  for  61 1 1-T4  aluminum. 
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Fig.  7.  Anisotropy  of  flow  stresses  (normalized  by  uniaxial  tension  about  the  x-axis)  and  r-values  for  6111-T4, 
measured  and  predicted  by  CPB06ex2  (k  =  O.fc'  =  0). 


Fig.  8.  Plane-stress  yield  surface  (oxy  =  0)  for  2090-T3  aluminum  sheet  according  to  the  CPB06ex2  yield  criterion 
( k  =  0 =  0)  yield  surface  and  experimental  data  (tension  only).  Data  after  Barlat  et  al.  (2005). 

(a  =  12).  Furthermore,  a  comparison  between  the  45°  yield  stress  values  predicted  using 
CPB06ex2  and  the  VPSC  values  reported  by  Barlat  et  al.  (2005)  is  given  in  Table  11.  Figs. 
8  and  9  illustrate  the  ability  of  the  model  to  represent  the  material  with  a  high  degree  of 
accuracy.  The  variation  in  r-values  with  6  peaks  between  0°  and  360°  indicates  that  this 
model  is  likely  to  be  suitable  for  the  prediction  of  six  or  more  ears  in  cups  drawn  from 
circular  blanks. 
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Fig.  9.  Anisotropy  of  flow  stresses  (normalized  by  uniaxial  yield  stress  along  the  rolling  direction)  and  /--values 
for  2090-T3,  measured  and  calculated  with  CPB06ex2  yield  criterion  (k  =  0,kl  =0).  Data  after  Barlat  et  al. 
(2005). 


Table  10 

CPB06ex2  yield  criterion  parameters  for  2090-T3  Al  alloy  (tension  only) 


c,  1 

C12 

Cn 

C22 

^23 

C33 

C44 

C55 

c66 

0.476 

-0.724 

-1.209 

-0.987 

-2.098 

-3.070 

0.563 

2.008 

1.117 

c' 

U11 

c' 

'-12 

c' 

c13 

c' 

*-22 

c' 

'-'23 

c' 

c33 

c' 

'-'44 

c' 

c55 

c' 

'-'66 

0.476 

-0.763 

1.143 

0.095 

-0.608 

0.715 

1.940 

1.032 

2.028 

Table  1 1 

Comparison  between  theoretical  normalized  yield  stress  values  for  2090-T3  Al  alloy 


Y-Z  45°  Tension 

Z-X  45°  Tension 

Y-Z  simple  shear 

Z-X  simple  shear 

VPSC 

0.90 

0.89 

0.47 

0.47 

CPB06ex2 

0.9 

0.89 

0.49 

0.47 

VPSC  results  after  Barlat  et  al.  (2005). 


The  last  example  pertains  to  Y350  Mpa  cold  rolled  high  strength  steel  (HSS).  The 
experimental  data  consisting  of  tensile  yield  stresses  and  r-values  along  seven  directions 
in  the  plane  of  the  sheet  (i.e.  from  rolling  to  transverse  directions  in  15°  increments)  were 
reported  in  Hu  (2005).  For  this  material  the  exponent  a  =  5.  The  values  of  the  coefficients 
of  anisotropy  coefficients  involved  in  the  CPB06ex2  yield  criterion  are  given  in  Table  12. 
Fig.  10  shows  the  plane-stress  yield  surface  calculated  with  CPB06ex2  using  the  experi¬ 
mental  yield  stresses.  Fig.  1 1  compares  the  experimental  and  predicted  anisotropic  prop¬ 
erties  of  the  steel  sample.  An  excellent  agreement  with  the  experimental  data  is  obtained. 
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Table  12 

CPB06ex2  yield  criterion  parameters  for  Y350  MPa  steel 


Cn 

Cl2 

Cn 

C22 

C23 

C33 

c66 

1.761 

0.630 

0.475 

1.712 

0.229 

1.664 

1.499 

c' 

'-'11 

c' 

c12 

c' 

'“'13 

c' 

'“'22 

c' 

^23 

c' 

u33 

c' 

c66 

1.761 

2.913 

2.661 

4.616 

3.358 

4.352 

0.761 

Fig.  10.  The  CPB06ex2  (k  =  0,12  =  0)  yield  surface  for  plane  stress  in  the  X-Y plane  for  F1SS  steel.  Experimental 
data  from  Hu  (2005). 


Fig.  11.  Anisotropy  of  flow  stresses  (normalized  by  uniaxial  tension  about  the  x-axis)  and  r-values  for  HSS  steel 
measured  and  predicted  by  CPB06ex2  (Ar  =  0).  Data  after  Hu  (2005). 


6.  Conclusions 

In  this  paper,  yield  functions  describing  the  anisotropic  behavior  of  textured  metals 
were  proposed.  These  yield  functions  are  extensions  to  orthotropy  of  the  isotropic  yield 
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function  proposed  by  Cazacu  et  al.  (2006).  Anisotropy  is  introduced  using  linear  transfor¬ 
mations  of  the  stress  deviator.  It  is  shown  that  the  anisotropic  yield  function  CPB06ex2 
which  involves  18  anisotropy  coefficients  (two  linear  transformations)  represents  with 
great  accuracy  both  the  tensile  and  compressive  anisotropy  in  yield  stresses  and  /-values 
of  materials  with  hep  crystal  structure  and  of  metal  sheets  with  cubic  crystal  structure. 
Furthermore,  it  was  demonstrated  that  it  can  describe  very  accurately  the  anisotropic 
behavior  of  metal  sheets  whose  tensile  and  compressive  stresses  are  equal. 

It  was  shown  that  the  accuracy  in  the  description  of  the  details  of  the  flow  stress  and  r- 
values  anisotropy  in  both  tension  and  compression  can  be  further  increased  if  more  than 
two  linear  transformations  are  included  in  the  formulation.  If  the  in-plane  anisotropy  of 
the  sheet  is  tension  and  compression  is  not  very  strong,  the  yield  criterion  CPB06ex2  pro¬ 
vides  a  very  good  description  of  the  main  trends. 

Appendix  A.  Derivatives  of  the  CPB06ex2  function 


In  the  current  work,  we  assume  an  associated  flow  rule.  Let  denote  by  r0  the  Lankford 
coefficient  (width  to  thickness  strain  ratios)  under  uniaxial  tensile  or  compressive  loading 
in  a  direction  at  angle  6  with  the  rolling  direction  in  the  plane  {xy).  Hence: 


sin2  0  eE  -  sin(20)  £ + cos2  6  £ 


re  = 


6 F 
da xx 


daX} 
8 F 
8(7yi/ 


(A.1) 


where  is  given  by  Eq.  (8).  The  derivatives  of  F  reduce  to  the  following  form  for 

uniaxial  tension  in  the  x-y  plane,  and  are  required  for  the  calculation  of  r0: 


8 F  8Xm  8Zkl  6 F  8I’m  81'kl 

nj  81  m  SI,.,  8a ij  8Z’m  8 1’kl  8a  tj 

The  non-zero  terms  in  (A. 2)  are  given  by 
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Relations  similar  to  (A. 3)  express  each  term  as  a  function  of  the  corresponding  primed 
variable. 
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